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SparkJet  Actuators  for  Flow  Control 


Executive  Summary 

The  generation  of  plasmadynamic  discharges  for  near-surface  aerodynamic  applications 
holds  the  promise  of  providing  rapid  and  flexible  steering  control  for  advanced  high-speed  flight 
vehicles.  Effective  manipulation  of  a  flow  field  can  lead  to  a  number  of  significant  benefits  to 
aerospace  vehicle  systems,  including  enhanced  performance,  maneuverability,  payload  and 
range,  as  well  as  lowered  overall  cost.  These  macro  benefits  are  directly  achievable  through  the 
application  of  flow  control  technology  to  impact  fluid  phenomena  such  as  transition,  turbulence, 
and  flow  separation  on  a  micro-scale.  Practical  application  of  active  flow  control  is  dependent 
upon  the  development  of  robust  actuators. 

Under  AFOSR  sponsorship  (AFOSR  Grant  Number  FA9550-04-1-0095),  The  Johns 
Hopkins  University  Applied  Physics  Laboratory  (JHU/APL)  is  investigating  a  promising  cavity 
device  for  high-speed  flow  control  called  the  SparkJet  actuator.  This  actuator,  which  produces  a 
synthetic  jet  with  high  exhaust  velocities,  holds  the  promise  of  manipulating  high-speed  flows 
without  active  mechanical  components.  Computational  and  experimental  techniques  are  being 
used  to  investigate  the  operating  characteristics  of  this  cavity  device.  Computational  Fluid 
Dynamics  (CFD)  modeling  is  used  to  investigate  the  fundamental  physics  and  flow  interactions 
within  and  near  the  device.  CFD  parametric  studies  characterized  the  performance  attributes  of 
the  mm-scale  device  as  a  function  of  orifice  size,  chamber  volume,  and  energy  deposited. 
Ongoing  experimental  efforts  to  assess  the  latest  actuator  design  include  the  application  of  high- 
resolution  particle  image  velocimetry  to  quantify  quiescent  air  operation,  and  the  use  of  a 
miniaturized  thrust  stand  to  measure  SparkJet  impulse  bit  data  and  to  determine  an  optimal  duty 
cycle.  Both  data  sets  are  being  used  to  calibrate  the  CFD  computations  and  to  determine  a 
confidence  level  in  the  prediction  capability. 

The  SparkJet  actuator  has  demonstrated  the  ability  to  produce  a  synthetic  jet  with  high 
exhaust  velocities  and  hence  holds  the  promise  of  manipulating  high-speed  flows  without 
moving  aerodynamic  structures.  Existing  piezo-mechanical  synthetic  jet  devices,  which  operate 
with  much  lower  exhaust  velocities,  have  been  shown  to  be  effective  for  low-speed  flow  control, 
but  they  cannot  develop  sufficient  driving  pressure  to  penetrate  supersonic  boundary  layers.  The 
envisioned  long-term  payoff  is  the  development  of  a  micro-actuated  SparkJet  flow  control 
system  for  enhanced  aerodynamic  performance  of  small,  highly  maneuverable,  high-speed 
vehicles. 

Personnel  Involved  (Duration  of  Grant) 

JHU/APL:  Bohdan  Z.  Cybyk  (aerospace  engineer),  Kenneth  R.  Grossman  (electrical  engineer), 
H.  Bruce  Land  III  (electrical  engineer),  Daniel  H.  Simon  (aerospace  engineer),  Jordan  Wilkerson 
(mechanical  engineer) 

JHU:  Jun  Chen,  post-doc,  Mechanical  Engineering  Department  (Oct  2004  through  Jul  2005) 
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1.  Introduction 

Active  flow  control  is  regarded  as  an  enabling  technology  for  many  advanced  air  vehicle 
concepts  under  consideration.  Effective  manipulation  of  a  flow  field  can  lead  to  a  number  of 
significant  benefits  for  aerospace  vehicles,  including  enhanced  performance,  maneuverability, 
payload,  and  range,  as  well  as  lowered  overall  cost.  These  macro  benefits  are  directly  achievable 
through  the  application  of  flow-control  technology  to  impact  fluid  phenomena  such  as  transition, 
turbulence,  and  flow  separation  on  the  micro  scale.1,2,3,4,5,6  Organizations  such  as  the  United 
States  Air  Force  and  NASA  continue  to  investigate  the  potential  advantages  of  active  flow 
control  over  more  traditional  aerodynamic  techniques.7,8  A  need  exists  for  a  robust  flow  control 
actuator  that  is  capable  of  penetrating  supersonic  boundary  layers. 

Several  classes  of  micro-actuators  are  under  investigation  for  flow  control  applications 
associated  with  aerospace  vehicle  systems.  Most  use  mechanical  deflection  of  control  surfaces, 
mass  injection,  or  synthetic  jets.  Others  manipulate  electromagnetic  fields  in  an  attempt  to 
control  flow.9, 10,11  Mechanical  actuators  include  electroactive  polymers,  shape  memory  alloys, 
electro-active  ceramics,  and  Micro-Electro-Mechanical  Systems  (MEMS)  flaps.12,13,1  These 
actuators  control  flow  by  movement  of  a  control  surface  or  manipulation  of  surface  texture. 
Mass  injection  devices  include  combustion-driven  jet  actuators,  which  burn  a  gaseous  fuel-air 
mixture.15  A  small  chamber  is  filled  with  a  combustible  mixture  and  then  ignited,  resulting  in 
high  pressures  inside  the  chamber  and  mass  expulsion  through  the  chamber  orifice.  Synthetic  jet 
actuators  are  fluidic  control  devices.16  They  typically  use  a  piezo-electric  diaphragm  in  a  cavity 
opposite  an  orifice.  The  oscillatory  motion  of  the  diaphragm  alternately  decreases  the  cavity 
volume,  expelling  gas,  then  increases  the  cavity  volume,  refilling  the  cavity  with  gas.  Synthetic 
jet  actuators  transfer  momentum  into  the  external  system  without  net  mass  transfer. 

Each  of  the  aforementioned  actuators  has  significant  limitations  for  application  to 
supersonic  flow  control.  Shape  memory  alloys  have  high  power  requirements,  poor  cycle  life, 
and  are  not  readily  scalable.  Electro-active  ceramics  produce  very  small  strain.  MEMS  devices 
are  highly  complex  and  expensive.  Mass  injection  devices  require  fuel  sources  and  mechanical 
valving.  Finally,  diaphragm-driven  synthetic  jets  produce  relatively  low  jet  velocities.  Under 
Air  Force  Office  of  Scientific  Research  (AFOSR)  sponsorship,  The  Johns  Hopkins  University 
Applied  Physics  Laboratory  (JHU/APL)  is  investigating  a  system  for  high-speed  flow  control 
called  the  SparkJet.  The  SparkJet,  which  produces  a  synthetic  jet  with  high  exhaust  velocities, 
holds  the  promise  of  manipulating  high-speed  flows  without  moving  structures.  The  overall 
objective  of  the  three-year  research  grant  is  to  develop,  assess,  and  mature  the  SparkJet  actuator 
concept  for  active  supersonic  flow  control.  Computational  and  experimental  techniques  are 
being  used  concurrently  to  quantify  the  performance  characteristics  of  the  pulsating  flow  from  a 
cavity  version  of  the  device.  Analyses  and  demonstrations  performed  to  date  indicate  that  the 
SparkJet  concept  has  the  potential  of  generating  exhaust  streams  that  can  penetrate  supersonic 
(as  well  as  subsonic)  boundary  layers  without  the  need  for  active  mechanical  components. 

2.  Nomenclature 

D  =  distance  between  collocated  actuators 

Dci,  =  chamber  diameter 
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=  chamber  orifice  diameter 
=  chamber  internal  energy,  Young’s  modulus 
=  resonant  frequency 
=  light  intensity 
=  impulse  bit 
=  spring  constant 
=  beam  length 
=  Mach  number 
=  mass 

=  effective  system  mass 
=  initial  effective  mass 
=  energy  deposited 
=  temperature 
=  time,  beam  thickness 
=  chamber  volume 
=  capacitor  potential 
=  velocity  in  the_p  direction 
=  beam  width 
=  amplitude,  distance 
=  chamber  axis 
=  flow  angle 

=  wavelength  of  the  laser  source 
=  effective  natural  frequency 
=  damped  effective  natural  frequency 
=  effective  damping  coefficient 


3.  Actuator  Description  and  Operation 

The  SparkJet  is  a  completely  solid-state  device  that  consists  of  a  small  chamber  with 
electrodes  and  a  discharge  orifice.17  High  chamber  pressure  is  generated  by  rapidly  heating  the 
gas  inside  the  SparkJet  using  an  electrical  discharge.  The  pressure  is  relieved  by  exhausting  the 
heated  air  through  an  orifice. 

A  single  cycle  of  SparkJet  operation  consists  of  three  distinct  stages:18  energy  deposition, 
discharge,  and  recovery  (Fig.  1).  The  cycle  is  repeated  continuously  to  create  a  synthetic  jet 
capable  of  causing  a  macro-scale  flow  effect 
when  positioned  in  a  region  of  high 
receptivity.  The  device  transfers  momentum 
with  the  external  surroundings  without  net 
mass  transfer.  The  amount  of  mass  expelled 
during  each  cycle  is  large  compared  to  that 
of  piezo-electric  devices. 

In  active  flow  control,  disturbances 
produced  by  such  a  cavity  device  are  used  to 
take  advantage  of  fluid  phenomena  such  as 


Stage  1:  Stage  2:  Stage  3: 

Energy  Discharge  Recovery 

Deposition 
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Figure  1.  SparkJet  cycle. 


SparkJet  Actuators  for  Flow  Control 
Final  Report 


Enclosure  to: 
GED-L-07-0137 


transition,  turbulence,  and  flow  separation.  Conversely,  these  disturbances  can  also  be  used  to 
study  the  development  of  such  phenomena.  For  example,  electrical  discharges  in  a  recessed 
cavity  have  been  used  in  experiments  to  generate  artificial  disturbances  for  the  study  of  laminar- 
to-turbulent  transition.19,20  The  disturbance  frequencies  were  at  levels  typically  needed  for  active 
flow  control. 

Figure  2  presents  a  cross  section  of  the  SparkJet  chamber.  The  device  consists  of  a 
ceramic  component  with  a  cathode  that  determines  the  chamber  volume.  The  ceramic 
component  is  inserted  between  a  stainless  steel  holder  and  a  threaded  brass  cap  with  a  jewel 
orifice  installed  in  the  face  (Fig.  3).  This  SparkJet  design  has  several  significant  advantages  over 
previous  designs:  it  is  more  robust,  allowing  for  reliable  testing;  it  provides  a  precise  volume  and 
orifice;  and  it  allows  for  greater  flexibility  in  testing  and  ease  of  use. 


Anode 


Figure  2.  Schematic  of  the  two-  Figure  3.  Individual  SparkJet  components 
electrode  SparkJet  actuator.  (inset  shows  assembled  unit). 


The  global  characteristics  of  a  single  discharge  cycle  have  been  determined  using 
analytical,  computational,  and  experimental  methods.  These  efforts  are  discussed  in  the  next 
three  sections. 

4.  First-Order  Model 

A  first-order  model  of  the  SparkJet  operation  has  been  developed  using  the  energy 
equation  applied  to  a  control  volume  enveloping  the  chamber: 


'qpdV  -  <£jV  dS-^PU-dS 
V  s  s 


(1) 
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In  this  equation  q  =  heat  added  per  unit  mass  to  the 
control  volume  and  q'=  heat  lost  from  the  control 
volume  through  its  surface.  (a) 

A  detailed  discussion  of  this  analytical 
model  can  be  found  in  Appendix  A  as  well  as  in 
Grossman  et  al.  (2003).  Using  this  analytical 
model,  basic  operating  characteristics  of  the 
SparkJet  device  can  be  discerned.  For  instance, 
many  aspects  of  the  operation  are  controlled  by  the 
ratio  of  the  energy  deposited  (Q)  to  the  initial 
internal  energy  (E)  of  the  gas  within  the  control 
volume.  As  seen  in  Fig.  4(a),  the  velocity  of  the 
discharge  jet  (Uex)  ranges  from  approximately  600 
to  950  m/s  at  the  beginning  of  the  fluid  discharge 
stage  for  sonic  orifice  flow.  The  amount  of  the  .  . 
initial  mass  expelled  and  the  fraction  of  deposited 
energy  expelled  are  shown  in  Fig.  4(b)  and  Fig. 

4(c),  respectively.  At  high-energy  deposition 
levels,  the  SparkJet  is  very  efficient  in  transferring 
electrical  energy  into  kinetic  energy  of  the  injected 
mass. 

5.  Computational  Predictions 

Both  2D  and  3D  time-dependent  Navier-Stokes  calculations  of  the  SparkJet  operation 
have  been  conducted  using  the  commercial  Computational  Fluid  Dynamics  flow  solver  CFD++. 
Figure  5  presents  a  1200  ps  time  history  of  the  temperature  and  velocity  fields  predicted  by  CFD 
within  and  above  the  SparkJet  chamber.  In  this  figure,  the  color  scale  progresses  from  blue 
(min)  to  red  (max)  corresponding  to  temperatures  from  300  to  1000  K  and  velocities  from  -100 
to  400  m/sec.  The  rectangular  zone  that  represents  the  annular  energy  deposition  region  in  3D  is 
visible  in  the  t=0  snapshot.  The  first  -1000  ps  period  represents  the  SparkJet  discharge  stage 
(stage  2)  of  device  operation,  followed  by  the  first  200  ps  of  the  recovery  stage  (stage  3).  Figure 
6  presents  the  time  histories  of  average  (area-weighted)  temperature  and  velocity  across  the 
orifice  exit  plane.  Note  that  substantial  fluctuations  exist  in  the  averaged  quantities  across  the 
orifice  exit  plane. 

A  parametric  CFD  study  was  conducted  to  quantify  the  sensitivities  of  single-pulse 
device  performance  on  orifice  diameter,  Dex,  chamber  volume,  Vch,  and  energy  deposited,  Q,  for 
the  SparkJet  firing  into  quiescent  air.  The  range  of  parameters  investigated  as  part  of  the  study 
included  orifice  diameters  of  5,  13,  and  26  mil  (0.127,  0.33,  0.66  mm),  chamber  volumes  of 
2.877E-8,  4.056E-8,  and  5.234E-8  m3,  and  energy  deposition  values  of  0.0126,  0.0302,  0.0478, 
and  0.0654  J.  The  effects  of  increasing  Q  on  total  impulse  are  shown  Fig.  7.  With  a  constant 
chamber  volume,  increasing  Q  increases  the  total  impulse.  Furthermore,  a  larger  diameter  orifice 
produces  a  larger  impulse,  indicating  that  viscous  effects  play  a  significant  role  for  the  range  of 
device  sizes  investigated.  Therefore,  a  SparkJet  can  be  calibrated  for  a  specific  application  by 


Q/E 


2.0  4.0  6.0  8.0 

Q/E 


Figure  4.  Effect  of  energy  deposition 
on  SparkJet  operating  characteristics 
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varying  the  energy  deposition  and  orifice  diameter.  CFD  simulations  of  fully-developed  Mach  3 
flow  over  a  surface-mounted  SparkJet  indicate  that  a  single  discharge  is  sufficient  to  penetrate  a 
supersonic  boundary  layer.21  SparkJet  operation  in  synthetic  jet  mode  (with  target  frequencies 
up  to  1  KHz)  significantly  increases  the  effect  of  the  effluent  on  the  surroundings. 
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Figure  5.  Temperature  and  velocity  time  histories  predicted  by  CFD  after 
single  energy  deposition  pulse  at  t=0. 
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Figure  6.  Average  temperature  and  velocity  time  histories  at  the  orifice  exit. 
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Figure  7.  Effect  of  energy  deposition  on  impulse. 


6.  Experimental  Results 

Time-sequenced  Schlieren  photographs,  shown  in  Fig. 
8,  provide  visual  verification  of  device  operation.  Although  it 
is  possible  to  obtain  some  quantitative  data  from  Schlieren 
photographs,  the  inherent  uncertainty  in  the  results  is 
relatively  high.  In  order  to  acquire  higher  accuracy  data, 
detailed  experimental  analyses  of  the  SparkJet  flow 
characteristics  are  being  explored  using  advanced  flow 
visualization  techniques.  Device  operation  is  being  examined 
using  high-resolution  particle  image  velocimetry  at  The 
Laboratory  for  Experimental  Fluid  Dynamics  at  the  Johns 
Hopkins  University.  Furthermore,  a  miniaturized  thrust  stand 
designed  to  support  the  development  of  a  micro-Pulsed 
Plasma  Thruster  concept  under  JHU/APL  IR&D  funding  is 
being  used  to  measure  SparkJet  impulse  bit  data  and  to 
determine  an  optimal  duty  cycle.  Both  data  sets  will  be  used 
to  calibrate  the  CFD  simulations  and  to  determine  a 
confidence  level  in  the  prediction  capability. 

6.1  Actuator  Redesign 

Initial  SparkJet  proof-of-principle  designs  and 
experiments  were  based  upon  the  simple  two-electrode  design 
shown  in  Figs.  2  and  3.  While  a  two-electrode  device  is 
relatively  simple  mechanically  and  electronically,  it  presents 
issues  operationally.  Operation  is  dependent  upon  the 
spontaneous  breakdown  of  the  air  for  creation  of  the  spark. 
The  breakdown  potential  will  vary  with  oxidation  of  the 
electrodes,  the  electrode  gap,  and  the  gas  pressure  within  the 
chamber.  Since  the  breakdown  is  spontaneous,  coordinating 
the  breakdown  with  diagnostic  equipment  is  difficult. 


Figure  8.  Schlieren  photo¬ 
graphs  of  a  SparkJet 
discharge. 
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Triggering  requirements  for  high-resolution  PIV 
studies  (discussed  in  the  next  section)  necessitated  a 
redesign.  The  design  has  now  transitioned  to  a  three- 
electrode  configuration  (Fig.  9)  to  reduce  the  size  of  the 
power  supply,  allow  easy  synchronization,  and  increase 
reproducibility  and  reliability.  The  new  devices  are 
machined  from  Macor,  a  machinable  glass  ceramic  from 
Corning.  Macor  was  chosen  because  of  its  insulative 
properties  and  the  ease  with  which  it  can  be  machined  with 
normal  tooling,  enabling  mass  production  as  cast  units. 

The  ability  to  control  the  timing  and 
synchronization  of  the  discharge  is  inherent  in  the  three- 
electrode  design,  allowing  experimental  measurements  to 
be  timed  precisely  with  device  operation.  The  three-electrode  device  has  a  ground  electrode,  a 
main  electrode,  and  a  trigger  electrode.  Device  discharge  occurs  only  when  the  trigger  voltage  is 
activated.  Thus,  the  voltage  and  energy  in  the  discharge  can  be  adjusted  over  a  wide  range 
without  causing  spontaneous  discharge  to  occur.  The  trigger  electrode  produces  a  kilovolt  pulse 
that  ionizes  the  air  gap,  allowing  the  energy  from  the  main  capacitor  to  discharge  between  the 
main  electrodes.  The  main  capacitor  needs  to  be  charged  only  to  hundreds  of  volts  rather  than 
the  kilovolts  needed  in  the  two-electrode  design,  resulting  in  significant  savings  in  power  supply 
cost,  weight,  and  volume. 

6.2  High-Resolution  PIV 

6.2.1  Setup  and  Instrumentation 

In  order  to  obtain  non-intrusive  field  measurements  needed  to  characterize  the  SparkJet 
discharge,  a  high-resolution  PIV  system  is  employed.22  The  SparkJet  device  is  mounted  inside 
and  on  the  top  center  of  an  acrylic  chamber  constructed  specifically  for  this  study.  The  test 
chamber,  shown  in  Figs.  10  and  1 1,  has  inner  dimensions  of  375(L)  x  250(W)  x  350(H)  mm.  It 

sits  on  an  optical  frame  with  PIV  hardware 
deployed  from  the  outside. 

The  test  chamber  is  seeded  with  smoke 
particles  (1-3  pm  in  diameter)  generated  by  a  Lite 
F/X  fog  machine  from  an  inlet  in  the  bottom 
corner  of  the  chamber.  A  vertical  baffle  plate  is 
set  up  inside  the  chamber  to  prevent  big  smoke 
particles  from  entering  the  sample  area.  A  small 
electrical  fan  with  DC  speed  control  introduces  a 
weak  circulating  flow  pattern  inside  the  chamber 
to  agitate  the  smoke  particles.  Only  small 
particles  enter  the  sample  area  while  large 
particles  are  stranded  near  the  bottom  and  far 
from  the  sample  area.  The  residual  flow  caused 


Figure  10.  Experimental  setup  formicro- 
PIV  studies  of  the  SparkJet  device. 


Figure  9.  Individual  SparkJet 
unit  (with  electronics  and  power 
supply). 
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by  the  agitation  process  is  very  small  compared  to  the  jet  velocity  and  thus  is  ignored.  The 
smoke  generator  ejects  smoke  particles  into  the  test  chamber  only  during  the  experiment 
preparation  stage.  Preliminary  tests  show  that  the  smoke  particles  remain  uniformly  suspended 
for  approximately  20  minutes.  The  internal  chamber  of  the  SparkJet  device  is  not  seeded  with 
particles. 

The  beam  from  the  PIV  light  source,  a  dual-head  Nd-YAG  laser  (532  nm,  120  mJ/pulse 
maximum  energy),  is  expanded  to  a  500  pm  thick  sheet  along  the  central  X1-X2  plane.  The  25  x 
25  mm  sample  area  is  located  directly  underneath  the  jet  exit.  The  images  are  recorded  using  a 
Kodak  ES-4.0  digital  camera,  with  a  resolution  of  2048  x  2048  pixels,  and  a  Nikon  105mm 
Nikkor  lens.  The  camera  is  mounted  on  a  3D  translation  stage  for  fine  position  adjustment.  The 
power  supply  of  the  SparkJet  device,  the  laser  and  camera  controllers,  and  the  data  acquisition 
system  are  all  synchronized  with  each  other,  as  shown  in  Fig.  12. 

The  camera  records  two  consecutively  exposed  images  on  two  separate  frames.  The  time 
interval  between  laser  pulses  ranges  from  2  to  10  ps,  depending  on  the  timing  phase  of  the  jet 
evolution,  and  is  too  short  to  be  controlled  accurately  by  only  setting  the  timing  of  firing  flash 
lamps.  Instead,  the  flash  lamps  of  the  two  lasers  are  fired  at  the  same  time,  and  the  delay 
between  the  two  laser  pulses  is  established  by  triggering  the  Q-switches  of  the  laser  heads  at 
different  times.23  Since  there  are  slight  variations  in  the  response  of  the  two  laser  heads  to  the 
trigger  signal,  the  exact  laser  timing  is  monitored  using  a  fast-response  photodiode  connected  to 
a  2.5  GS/s  oscilloscope.  The  trigger  signals  for  the  flash  lamps  and  Q-switches  are  generated  by 
a  4-channel  Stanford  Research  Delay  Generator.  This  setup  produces  a  timing  sequence  of  laser 
pulses  with  an  accuracy  of  20  nanoseconds. 


100mm 


Smoke  Generator 


375mm 


Figure  1 1 .  Schematic  of  the  SparkJet  test  chamber.  Front  view  (left),  side  view  (right). 
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Figure  12.  Schematic  diagram  of  the  PIV  instrumentation  system. 


6.2.2  Data  Processing 

The  edge  portion  of  the  acquired  images  is  contaminated  by  surface  reflection  and  is 
discarded  before  data  processing.  The  data  processing  consists  of  image  enhancement,  followed 
by  cross-correlation  analysis  to  determine  the  velocity  distributions.24,25,26  The  size  of  the 
interrogation  window  is  48  *  48  pixels.  With  50%  overlap,  the  vector  spacing  is  300  pm.  A 
total  of  63  x  75  vectors  are  obtained  from  each  image  pair. 

6.2.3  Observations  and  Results 

The  SparkJet  device  being  investigated  has  an  orifice  diameter  of  Dex  =  0.33  mm  (0.013 
inches)  and  a  chamber  volume  Vc/,  =  4.228E-8  m3.  The  capacitor  potential  to  initiate  the  spark 
discharge  is  V,-  =  1000  volts.  The  SparkJet  is  operated  at  0.2  Hz  to  examine  the  characteristics  of 
a  single  pulse.  Figure  13  presents  typical  PIV  images  from  two  different  SparkJet  discharges 
under  the  same  experimental  conditions.  The  flow  structure  is  well  visualized  because  the 
ambient  air  within  the  test  chamber  contains  enough  smoke  particles  to  be  seen  in  the  recorded 
images,  while  the  jet  core  is  dark  due  to  insufficient  seeding.  The  jet  structure  exhibits  a  good 
repeatability,  making  it  feasible  to  study  the  evolution  of  the  jet  through  repeated  measurements 
of  different  SparkJet  discharges. 

Figure  14  shows  the  resolved  velocity  fields  at  different  snapshots  in  time  as  well  as  the 
vorticity  components  calculated  using  finite  differencing.  All  these  results  are  considered 
instantaneous  due  to  the  very  short  laser  pulse  (on  the  order  of  10  ns)  and  no  ensemble  averaging 
is  applied.  At  t  -  0,  the  spark  is  ready  to  discharge;  no  flow  motion  is  detected.  After  the  spark 
discharges,  a  puff  of  hot  plasma  is  ejected  into  the  PIV  test  chamber  through  the  orifice.  This 
puff  propagates  into  the  test  chamber  and  expands  in  all  lateral  directions,  introducing  a  strong 
leading  vortex  ring.  Behind  this  puff,  the  ambient  air  is  entrained  into  the  wake,  producing  a 
weak  trailing  vortex  ring.  Both  vortex  ring  structures  can  be  seen  clearly  at  t  =  25  ps.  Visible  in 
each  velocity  field  snapshot  is  the  wave  front  propagation  at  t  =  25,  50,  75,  and  150  ps  when  the 
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air  at  the  puff  interface  is  significantly  compressed.  At  t  =  ~200  ps  (not  shown),  the  measured 
velocity  of  the  entrained  fluid  reaches  a  peak  value  of  100  m/s.  The  actual  speed  at  the  jet  core  is 
believed  to  be  much  higher  due  to  insufficient  seeding  in  the  core  area  and  the  spatial-averaging 
effect  of  the  PIV  velocity  measurements.  As  the  mixing  between  the  jet  and  the  ambient  air 
continues,  increasingly  complex  vorticity  patterns  are  produced.  At  t-  500  ps,  the  jet  penetrates 
the  entire  depth  of  the  field  of  view  (25  mm).  The  jet  continues  to  propagate  further  into  the  test 
chamber  as  the  velocity  magnitudes  decay. 


Figure  13.  Particle  images  of  separate  SparkJet  discharges  (t  =  200  ps). 

6.2.4  Lessons  Learned 

Characterization  of  the  unique  SparkJet  device  using  PIV  has  proven  to  be  challenging. 
The  primary  hurdle  is  the  proper  seeding  of  the  internal  SparkJet  chamber.  In  most  previous  PIV 
investigations  of  jet  flows,  7  both  the  emanating  jet  and  the  ambient  fluid  are  seeded 
appropriately.  In  the  present  study,  only  the  ambient  fluid  has  been  seeded  because  the  chamber 
inside  the  SparkJet  device  is  very  small  and  a  tiny  orifice  is  the  only  access  point.  Consequently, 
the  actual  high-speed  jet  core  flow  has  yet  to  be  measured.  In  an  attempt  to  overcome  this 
problem  without  modifying  the  geometry  of  device,  we  are  exploring  the  possibility  of  seeding 
the  spark  chamber  directly  by  injecting  a  very  small  amount  (0.005ml)  of  methanol  solution  with 
SiC>2  particles  (mean  diameter  850  nm)  through  the  orifice.  The  volatilizing  methanol  is 
expected  to  leave  only  SiC>2  particles  in  the  spark  chamber.  Once  the  seeding  problem  is  solved, 
PIV  measurements  will  be  repeated  as  needed  for  statistical  analysis. 

Strong  electromagnetic  interference  (EMI)  from  the  SparkJet  discharge  had  been  an  issue 
with  the  sensitive  control  apparatus  for  the  Nd-YAG  laser  and  camera  control.  False  triggers 
were  common  and  had  prevented  collection  of  quantifiable  velocity  data.  The  SparkJet  power 
supply  redesign  (described  previously)  significantly  reduced  the  EMI,  as  did  shielding  of  the 
power  supply  enclosure  using  aluminum  foil. 

Glare  emitted  by  the  spark  was  also  of  concern  because  it  may  degrade  the  particle 
images.  However,  preliminary  tests  indicated  that  this  undesired  light  is  much  weaker  than  the 
laser  sheet.  Since  the  adverse  influence  on  particle  images  is  minimal,  an  optical  bandpass  filter 
initially  added  to  block  the  glare  is  no  longer  necessary. 
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Figure  14.  Time  history  of  the  resolved  PIV  velocity  fields  (left)  and  the  out-of-plane 
component  of  vorticity  (right). 

Limitations  of  the  available  double-exposure  camera  have  also  impacted  the  application 
of  PIV  to  the  Sparklet  actuator.  The  Kodak  ES-4.0  camera  is  capable  of  taking  PIV  images  with 
a  minimum  laser  pulse  delay  of  1  ps  for  cross-correlation  analysis.  At  present,  the  JHU 
Laboratory  of  Experimental  Fluid  Dynamics  is  purchasing  a  new  state-of-the-art  PIV  system 
with  digital  recording  capability  (minimum  laser  pulse  delay  200  ns)  that  is  expected  to  meet 
Sparklet  analysis  requirements. 
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6.3  Miniaturized  Thrust  Stand 

A  miniaturized  thrust  stand  is  under  development  at  JHU/APL  to  assist  with  the 
characterization  of  various  micropropulsion  devices.  The  stand,  pictured  in  Fig.  15,  consists  of  a 
small  3  x  3  cm  platform  suspended  by  four  rigid  cantilever  beams.  The  platform  and  beams  are 
fabricated  from  a  single  sheet  of  stainless  steel  shim  stock.  This  metal  structure  is  oriented 
vertically  and  clamped  at  the  top  and  bottom  in  a  Lexan®  fixture.  The  fixture  is  mounted  on  top 
of  a  custom-built  miniature  shaker  table  that  is  used  to  calibrate  the  thrust  stand.  The  shaker 
table  is  driven  by  a  single  piezoelectric  linear  actuator  from  ThorLabs.  A  SparkJet  actuator  can 
be  bonded  directly  to  this  movable  3  x  3  cm  platform.  Electrical  power  is  carried  to  the  device 
by  several  strands  of  thin  magnet  wire  (small  gauge  wire  with  an  electrically  insulating  shellac 
coating).  Impulse  imparted  from  the  actuator  causes  the  platform  to  vibrate  ever  so  slightly. 


Figure  15.  JHU/APL  miniature  thrust  stand,  (a)  Front,  (b)  Back. 


To  measure  the  amplitude  of  these  vibrations  a  Michaelson  interferometer  setup  is 
employed.  As  illustrated  in  Fig.  16,  the  beam  from  a  4.5  mW,  670  nm  laser  diode  module  from 
ThorLabs  is  passed  through  a  non-polarizing  beam  splitter.  One  portion  of  the  beam  then  passes 
through  a  quartz  window  into  the  vacuum  chamber,  reflects  off  a  mirror  mounted  on  the  thrust 
stand’s  moveable  platform,  and  returns  to  the  beam  splitter.  The  other  portion  of  the  beam  is 
reflected  back  to  the  beam  splitter  through  a  stationary  mirror.  The  beam  splitter  then  directs 
both  beams  towards  a  photodiode  detector  where  they  generate  an  interference  pattern.  As  the 
position  of  the  thrust  stand  platform  moves,  the  intensity  of  the  light  at  the  detector  varies  as 


7(0  = 


'^r1+cJ*Ep>j 


(2) 


where  x(t)  is  the  position  of  the  mirror  and  A  is  the  wavelength  of  the  laser  source.  Note  that  the 
output  is  sinusoidal  and  repeats  every  time  the  platform  moves  by  a  14  of  a  wavelength.  To 
minimize  background  vibrations  due  to  environmental  effects,  the  vacuum  chamber  and 
interferometer  have  been  assembled  on  an  optical  bench  with  pneumatic  isolation  legs. 


16 


Sparklet  Actuators  for  Flow  Control 
Final  Report 


Enclosure  to: 
GED-L-07-0137 


Floating  Optical  Bench 


Test 

Chamber  Window 


SparkJet 


C - 

J 

$ 


Stationary  Mirror 
Beam  Splitter 

7  n 

Hj  Detector 


Footprint  Platform 


Figure  16.  Laser  interferometer  arrangement  for  measuring  thrust  stand  vibrations. 


In  addition  to  a  high  natural  frequency,  the  rigid  design  of  the  steel  platform  and  affixed 
beam  structure  is  intended  to  have  several  other  effects.  First,  the  impulse  imparted  by  a  test 
device  should  cause  a  small  deflection  in  each  of  the  beams  such  that  they  behave  like  simple 
linear  springs.  Second,  the  stiffness  of  the  beams  should  overwhelm  that  of  the  power  supply 
wires  and  make  the  system  simple  to  model.  And  third,  the  beams  can  limit  the  deflection  of  the 
platform  to  less  than  !4  of  the  interferometer  laser  wavelength  (to  avoid  multiple  fringe  patterns 

in  the  interferometer  data).  Larger  deflections  are  tolerable  but  call  for  a  more  sophisticated  data 

28 

reduction  process,  such  as  that  employed  by  Cubbin. 

Since  the  thrust  stand  has  been  designed  to  behave  like  a  damped  spring-mass  system,  its 
response  to  an  applied  impulse  is  given  by 


x(t)  =  e  — Lkii — sjnfftfy  efft) 


(3) 


where  ha  is  the  impulse  bit,  me/f  is  the  effective  system  mass,  %eff  is  the  effective  damping 
coefficient,  ft^is  the  effective  natural  frequency,  and  a>dejr  =  is  the  damped  effective 

natural  frequency.  If  the  damping  coefficient  is  small,  the  system  will  essentially  respond  to  an 
impulse  bit  with  a  regular  sinusoidal  response.  The  impulse  bit  can  then  be  inferred  by 
examining  the  amplitude  of  these  oscillations  using  the  approximation 

hu  x-meff-  cod  eff  (4) 


where  Ax  is  the  amplitude  of  the  thrust  stand  platform  vibrations.  The  only  other  system 
information  required  is  the  effective  mass  and  damped  natural  frequency. 

Each  of  the  beams  attached  to  the  moveable  platform  behaves  like  a  guided  cantilever 
whose  spring  constant  can  be  calculated  according  to  the  expression 


17 


SparkJet  Actuators  for  Flow  Control 
Final  Report 


Enclosure  to: 
GED-L-07-0137 


where  E  is  the  material  Young’s  modulus,  w  is  the  beam  width,  t  is  the  beam  thickness,  and  L  is 
the  beam  length.  This  system  has  four  beams  that  effectively  act  like  springs  in  parallel; 
therefore,  the  total  spring  constant  of  the  system  should  be  four  times  that  of  a  single  beam.  The 
effective  mass  of  the  system  should  be  close  to  the  actual  mass  of  the  vibrating  components 
(moveable  platform,  actuator,  main  capacitor,  and  mirror).  The  first  platform  and  beam  structure 
for  the  thrust  stand  was  designed  for  an  effective  mass  of  approximately  25  g.  The  beam 
dimensions,  w  =  0.25”  (6.35  mm),  t  =  0.008”  (0.20  mm),  and  L  =  0.33”  (8.3  mm),  were  then 
designed  to  provide  a  total  system  spring  constant  of  78,270  N/m  and  a  predicted  natural 
frequency  of  280  Hz.  Moreover,  this  design  was  expected  to  produce  vibrations  with  amplitudes 
ranging  from  40  to  175  nm  (a  comfortable  range  to  detect  with  an  interferometer  employing  a 
670  nm  red  laser  diode)  in  response  to  impulse  bits  ranging  from  2  to  7.5  p.N-sec. 

In  order  to  accurately  interpret  the  measurements  recorded  from  this  thrust  stand,  a 
suitable  calibration  process  must  be  applied.  The  conventional  approach  to  accomplish  this 
would  involve  the  application  of  a  series  of  known  loads  to  the  platform  and  measurement  of  the 
associated  deflections.  In  lieu  of  this  procedure,  an  unconventional  calibration  process  was 
developed  for  this  system.  As  mentioned  above,  the  thrust  stand  is  mounted  on  a  shaker  table 
that  can  impart  small  vibrations  controlled  by  a  standard  function  generator.  The  resonant 
frequency  of  the  platform/beam  structure  can  be  detected  by  varying  the  frequency  of  these 
vibrations  and  monitoring  the  amplitude  of  the  platform’s  oscillations.  Recall  from  Eq.  4  that 
this  quantity  alone  is  an  important  system  characteristic.  If  a  known  amount  of  mass,  Am,  is  then 
added  to  the  platform,  the  resonant  frequency  should  drop  in  a  predictable  manner.  Indeed  the 
new  resonant  frequency  of  the  system  should  behave  as 


In' 


m0,eff+Am 


(6) 


where  mo,eff  is  the  initial  effective  mass  (prior  to  the  addition  of  any  ballast).  This  expression  can 
be  rearranged  as 


A  m  =  - 


4  nl 


fd,eff  ) 


"0  ,eff 


(7) 


If  a  series  of  additional  weights  are  added  to  the  platform,  Eq.  7  tells  us  that  the  total  amount  of 
added  weight  should  vary  as  the  inverse  of  the  damped  natural  frequency  squared.  Therefore, 
these  two  quantities  can  be  plotted  and  fit  with  a  straight  line.  They-intercept  of  this  line  reveals 
the  effective  mass  of  the  platform/beam.  The  slope  of  this  line  also  reveals  the  system  spring 
constant.  The  spring  constant  is  not  explicitly  required  to  interpret  the  response  of  the  thrust 
stand  to  an  impulse  excitation;  however,  it  is  needed  to  interpret  the  steady  platform  deflection 
produced  in  response  to  a  constant  thrust. 

This  miniature  thrust  stand  was  used  to  gather  preliminary  impulse  bit  data  for  a  single 
SparkJet  device.  The  SparkJet  was  powered  by  a  220  nF  capacitor  rated  at  1000  V.  Impulse  bit 
data  was  collected  at  five  discrete  voltage  levels  ranging  from  520  to  900  V.  These  points 
correspond  to  energy  deposition  levels  of  30,  40,  54,  70,  and  89  mJ.  A  sample  of  this  data  is 
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shown  in  Fig.  17.  The  impulse  bit  data  appears  to  vary  linearly  with  the  deposited  energy.  The 
impulse  bit  data  exhibits  some  scatter  at  each  energy  level.  This  could  be  due  to  inherent 
variation  in  the  SparkJet  performance  or  the  measurement  process  itself.  Most  encouraging,  Fig. 
18  shows  virtually  no  difference  in  the  thrust  stand  response  between  1  Hz  and  20  Hz  SparkJet 
operation.  Additional  development  work  on  the  thrust  stand  is  proceeding  with  the  aim  of 
reducing  its  uncertainty. 


Figure  17.  Preliminary  SparkJet  impulse  bit  measurements. 


Time  (s) 


Figure  18.  Miniature  thrust  stand  response  during  SparkJet  operation. 


7.  Distributed  Actuator  Array  Model  Development 

Practical  application  of  the  SparkJet  technology  will  require  the  development  of 
distributed  actuator  arrays.  The  diminutive  size  (120-600  pm)  and  high  discharge  velocities 
(400-500  m/s)  of  the  SparkJet  actuator  make  it  very  difficult  to  evaluate  demonstration  units 
using  non-intrusive  measurement  techniques.  Hence,  modeling  must  be  used  to  evaluate  the 
effectiveness  of  distributed  actuator  arrays  and  their  potential  high-payoff  applications.  We  have 
begun  development  of  a  new  CFD-based  computational  methodology  that  will  enable  the 
accurate  and  cost-effective  simulation  of  problems  that  involve  distributed  actuator  arrays. 
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From  a  mathematical  standpoint,  CFD  formulations  will  rigorously  capture  the  flow 
control  effects  of  micro-scale  perturbations  on  the  surrounding  macro-scale  flow,  and  vice  versa. 
In  reality,  however,  these  spatially-challenged  solutions  require  extensive  computational 
resources  (both  in  solution  time  and  computer  memory)  that  render  such  an  approach  impractical 
for  more  than  a  few  perturbation  sources  operating  in  tandem.  This  is  because  modeling  the 
fluid  dynamics  within  a  micro-scale  perturbation  source  (e.g.  modeling  the  effects  of  the  energy 
deposition  within  the  micro-actuator  during  a  single  operating  cycle)  substantially  increases  the 
grid  and  computational  requirements  over  what  is  required  to  resolve  the  effects  of  the  imparted 
flow  disturbance  on  the  surrounding  flow.  Once  a  detailed  numerical  simulation  of  the 
perturbation  source  exists,  a  fluid  dynamic  model  of  the  operating  cycle  of  the  disturbance 
source  (a  virtual  source)  can  potentially  be  developed  and  implemented  via  boundary  conditions. 
This  is  precisely  the  approach  that  has  been  taken  by  others  evaluating  the  effectiveness  of 
synthetic  jet  actuators.29 

Like  the  synthetic  jet  actuator,  the  Sparklet  is  also  a  zero  net-mass-flux  device. 
However,  the  SparkJet  cycle  adds  an  additional  complication  that  invalidates  the  current 
boundary  condition  implementation  of  virtual  disturbance  sources.  From  a  modeling  standpoint, 
the  problem  is  that  the  SparkJet  recovery  stage  is  dependent  upon  the  external  flow  and  may  vary 
from  pulse  to  pulse.  Hence,  any  use  of  virtual  SparkJet  sources  should  capture  the  effect  of  the 
local  boundary  layer  flow  on  the  device  operation.  In  other  words,  the  virtual  sources  must  “self- 
adapt”  to  the  local  time-varying  environment. 

A  new  parametric  study  is  underway  to  evaluate  the  performance  characteristics  of 
collocated  pulsed  SparkJet  actuators  operating  in  a  Mach  3  flow  (Fig.  9).  The  purpose  of  this 
parametric  study  is  to  evaluate  the  effects  of  the  distance  between  the  actuators,  D,  and  incoming 
flow  angle,  a,  on  the  resultant  interacting  jets.  SparkJet  geometric  characteristics  (Fig.  10)  and 
the  amount  of  energy  deposited  are  being  held  fixed.  This  study  will  provide  guidelines  to  the 
development  of  self-adapting  virtual  SparkJet  sources. 


Figure  9.  Parametric  study  of  collocated  actuators.  Figure  10.  Collocated  actuator  grids. 
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8.  Summary 

Effective  manipulation  of  a  flow  field  can  lead  to  a  number  of  significant  benefits  to 
aerospace  vehicle  systems,  including  enhanced  performance,  maneuverability,  payload  and 
range,  as  well  as  lowered  overall  cost.  These  macro  benefits  are  directly  achievable  through  the 
application  of  flow-control  technology  on  the  micro  scale  to  impact  fluid  phenomena  such  as 
transition,  turbulence,  and  flow  separation.  Practical  application  of  active  flow  control  is 
dependent  upon  the  development  of  robust  actuators  that  are  reliable,  low  cost,  responsive,  and 
consume  little  power.  An  example  of  promising  actuator  technology  for  supersonic/hypersonic 
flow  applications  is  the  SparkJet  actuator  under  development  at  JHU/APL.  The  SparkJet  has 
demonstrated  the  ability  to  produce  a  synthetic  jet  with  high  exhaust  velocities,  and  hence  holds 
the  promise  of  manipulating  supersonic/hypersonic  flows  without  moving  aerodynamic 
structures.  Practical  application  of  the  SparkJet  technology  will  require  the  development  and 
demonstration  of  distributed  actuator  arrays  capable  of  operating  in  high-temperature 
environments. 

This  paper  provides  a  synopsis  of  SparkJet  technology  development  under  AFOSR  Grant 
Number  FA9550-04- 1-0095  and  introduces  an  approach  to  the  development  of  a  computational 
model  for  the  accurate  and  cost-effective  simulation  of  a  distributed  SparkJet  array  operating  in  a 
macro-scale,  high-speed  environment.  Combined  computational  and  experimental  techniques 
are  being  applied  to  investigate  the  operating  characteristics  of  a  cavity  SparkJet  device  with  the 
ultimate  goal  of  developing  a  micro-actuated  flow  control  system  for  enhanced  aerodynamic 
performance  of  high-speed  flight  vehicles. 
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Appendix  A:  First-Order  Sparklet  Model 


A  first  order  model  of  the  Sparklet  operation  has  been  developed  using  the  energy 
equation  applied  to  the  control  volume  shown  in  Fig.  Al: 


%qpdV-§q'dS-§PV-dS  =  $A 

1/  C  C  1/ 


(A.l) 


where  q  is  the  heat  added  per  unit  mass  to  the  control  volume  and  q'  is  the  heat  lost  from  the 
control  volume  through  the  surface  of  the  control  volume. 

m 


Figure  A.  1 .  Spark  Jet  configuration. 


Making  the  assumption  of  one-dimensional  flow  within  the  control  volume,  Eqn.  (A.l) 
simplifies  to: 

UpeV)  =  Q-QAsmrxE-mh,  (A.2) 

at 

where  p  -  density,  e  -  internal  energy,  Q  -  rate  of  energy  added  to  the  control  volume,  Q'  -  rate 
of  energy  lost  through  the  control  surface  (per  unit  area),  As  -  internal  surface  area  m  -  mass 
flow  exiting  control  volume,  and  ht  -  total  enthalpy  of  flow  exiting  control  volume. 

Stage  1:  Energy  Deposition 

Within  Stage  1,  an  electrical  discharge  is  created  within  the  Sparklet  volume.  The  energy 
addition  is  assumed  to  occur  fast  with  respect  to  both  heat  and  mass  loss,  so  Eqn.  (A.2)  becomes: 

PeV)=Q  (A. 3) 

at 

Assuming  that  the  energy  addition  is  constant  in  time,  the  fluid  can  be  represented  as  a 
calorically-perfect  gas  and  states  0  and  1  represent  the  conditions  at  the  beginning  and  end  of  the 
first  stage,  Eqn.  (A. 3)  can  be  written: 
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T,=T0+-O-  (A.4) 

CVM 

where  Mis  the  total  mass  inside  the  control  volume  and  Q  is  the  total  energy  added  to  the  control 
volume.  Since  the  gas  density  is  constant  during  this  stage,  the  pressure  at  the  end  of  Stage  1  can 
be  determined  from: 


P^PoR 


T0  + 


Q 


CVM 


(A. 5) 


Hence,  one  can  see  that  both  the  pressure  and  temperature  of  the  gas  within  the  volume  are 
increased  during  Stage  1 . 

Stage  2:  Fluid  Discharge 

Within  Stage  2,  the  hot,  high-pressure  gas  within  the  control  volume  discharges  through 
the  orifice  until  the  point  where  the  external  and  internal  pressures  are  equalized.  During  this 
stage,  no  energy  is  added  to  the  control  volume  and  heat  loss  through  the  wall  can  be  neglected, 
so  Eqn.  (A. 2)  can  be  written: 


—  (peV)  =  -mh,  (A.6) 

dt 

While  the  discharge  flow  is  choked  at  the  orifice  (i.e.,  while  P>\.89*P0),  the  mass  flow  can  be 
related  to  the  chamber  pressure  and  temperature  as  follows: 


(A. 7) 


m  = 


P  • 
—  m 

J 


PA,, 


where 


yP, 


yf  r+r 


y+ 1  ' 

2(r  i). 


If  state  2  is  the  point  in  the  expansion  process  where  Pi=Pq  and  the  expansion  is  assumed  to 
occur  isentropically,  the  time  required  for  the  fluid  discharge  can  be  determined  from: 


A  -7,  =■ 


f  P  °  A 
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T  =  T 

i2 


y± 

r 


(A.9) 


In  actual  operation,  the  time  required  for  Stage  2  will  be  somewhat  longer  than  that  estimated 
using  Eqn.  (A. 8)  since  the  flow  will  continue  to  discharge  past  the  point  where  the  orifice 
unchokes.  At  the  conclusion  of  Stage  2,  the  fluid  within  the  control  volume  is  at  ambient 
pressure,  but  with  a  temperature  that  is  above  ambient. 

Stage  3:  Recovery 

During  the  recovery  stage,  heat  is  lost  though  the  wall  of  the  Sparkjet.  As  the  heat  is  lost, 
fresh  air  is  drawn  into  the  control  volume  such  that  nearly  ambient  pressure  is  maintained.  The 
amount  of  energy  that  needs  to  be  transferred  to  the  walls  can  be  determined  from: 

A  £  =  p2V{e2-ei)  (A.  10) 


Assuming  that  the  heat  transfer  can  be  modeled  in  terms  of  a  film  coefficient: 

Q' 


h  = 


A  SURFACE  i^W  T) 


and  assuming  that  a  lumped  heat  capacity  system  model  can  be  used,  the  temperature  within  the 
control  volume  can  be  determined  from: 


T-T0 

t2-t0 


=  e 


-Ur 


(A.ll) 


where  T  = 


CvpV 

fa  A 

rtyi  SURFACE 


Thus  one  can  see  that  the  recovery  time  can  be  decreased  (hence  the  operating  frequency  can  be 
increased)  by  increasing  the  internal  surface  area  of  the  film  coefficient. 
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